Constraining vacuum gap models of pulsar radio emission using 

the intensity modulation index 

Janusz Gil^ & Fredrick A. Jenet^ 
ABSTRACT 

Recent observations suggest that the level of pulse-to-pulse intensity mod- 
ulation observed in a given radio pulsar may depend on its period and period 
derivative. Such a "modulation index relationship" (MIR) may be an important 
tool for determining the physical processes behind the radio emission. In the 
context of sparking gap models, the exact functional form of the MIR depends 
on the physical processes occurring on the surface of the neutron star in the 
region known as the "vacuum gap." Several possible vacuum gap models are 
studied here in order to determine the expected MIR for a given model. Current 
observations are consistent with two of the four models studied: the curvature 
radiation driven vacuum gap and the curvature radiation driven near threshold 
vacuum gap (CR-NTVG). It is shown that the inverse Compton scattering driven 
vacuum gap models are not supported by the current data. Given that the cur- 
rent data best supports the CR-NTVG model, its possible that all pulsars have 
strong (Ri 10^^ G) surface magnetic fields. 

Subject headings: pulsars: general 



1. INTRODUCTION 

The physical mechanism of the coherent pulsar radio emission has remained elusive since 
their discovery over 30 years ago. The observed high brightness temperatures together with 
enormous amount of phenomenology exhibited makes these sources very difficult to under- 
stand. However, it is generally accepted that pulsar radio emission is generated within a dense 
electron-positron plasma, created near the polar cap and flowing along open magnetic field 
lines. The observed pulse-to-pulse intensity modulation can arise from the time-dependent 
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lateral structure of this flow, probed once per pulsar period by the observer's line of sight. 
Ruderman & Sutherland (1975, RS75 henceforth) proposed a pulsar model in which the 
lateral structure was in the form of localized spark discharges in a charge depleted region 
near the stellar surface know as the "vacuum gap." Gil & Sendyk (2000, GSOO henceforth) 
explored the RS75 model in an attempt to relate the observed radio emission properties to 
a pulsar's period P and its derivative P. They postulated that the polar cap is populated 
as densely as possible with a number of sparks, each having a characteristic size as well 
as separation from adjacent sparks approximately equal to the height h of the vacuum gap 
acceleration region. This leads directly to the so-called complexity parameter 

equal to the ratio of the polar cap radius to the characteristic spark dimension. One can show 
that a is approximately equal to the maximum number of sparks across the polRx cap and 
thus the maximum number of sub-pulses that can appear in the pulse window. Therefore, a 
describes the complexity of single pulses and/or mean profile (see GSOO for details). 

The pulse-to-pulse intensity modulation can be quantified by the phase resolved modu- 
lation index defined as 

"^^'^^ ^ TJM^ ' 

where /s(0) is the pulsar intensity at pulse phase and the angle brackets represent averaging 
over the ensemble of single pulses (for details see Jenet, Anderson & Prince 2001). Within 
the framework of sparking gap model the observed pulse-to-pulse intensity modulation is 
due to the presence of a number sparks moving over the polar cap either erratically or 
in an organized manner (drifting). As the number of sparks increase, one expects to see 
less pulse-to-pulse intensity modulation. Hence, the modulation index (Eq. [1]) should be 
anticorrelated to the complexity parameter (Eq. [2]), to the extent that at very high values 
of a there should be no detectable intensity modulation (GSOO). 

In general, the complexity parameter (Eq. [1]) can be expressed in the form a oc 
Since the modulation index may be directly related to some unknown function of a, it is useful 
to use the Spearman rank ordered correlation (SROC) technique in order to determine if a 
correlation exists between a and m. The SROC technique is independent, up to a sign, of 
any monotonic function applied to m and/or a. In this case, only the ratio « = /3/ 7 is 
important. Jenet & Gil (2003a, JG03 henceforth) used the SROC technique to search for 
the proposed relationship, hereafter referred to as a modulation index relationship (MIR), 
using a sample of 12 pulsars. This preliminary investigation found evidence for a MIR with a 
between [—5, —2]. The most significant value of a = —2.7. In general, a MIR will provide an 
important constraint on any model of pulsar radio emission. Even its non-existence may be 
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an important piece of information. For the purposes of this letter, only spark gap models will 
be considered. These models predict the existence of a MIR together with a value of a. The 
specific value of a depends on the vacuum gap model in question. Thus, an observationally 
determined range of a. values will rule out various vacuum gap models. 

In the next section, four models of the vacuum gap acceleration region are studied in 
order to determine the expected value of a for each model. This letter is summarized in §3. 

2. MODELS OF VACUUM GAP ACCELERATION REGIONS 

All vacuum gap (VG) models assume that the free outflow of electrons or ions from the 
polar cap surface is strongly inhibited, which leads to formation of a charge depleted region 
just above the polar cap. These regions are conventionally called the vacuum gap even if the 
charge density is not exactly zero within it (for review see Gil, Melikidze & Geppert 2003). 
The charge depletion with respect to the corotational value (Goldreich & Julian 1969) results 
in an extremely high potential drop across the VG, which is discharged by photon induced 
pair creation in the strong and curved surface magnetic field. The details of this discharge 
depend on the strength and curvature of the surface magnetic field and on the physical 
process that generate the electron-positron pair plasma. 

The analysis presented here will calculate the complexity parameter for each VG model 
in question. In order to do this, both the vacuum gap height, /i, and the polar cap radius, r^, 
need to be estimated. Both of these quantities will depend on the structure of the magnetic 
field near the surface of the neutron star. 

Classically, the magnetic field is assumed to be dipolar. There is a growing evidence that 
the actual surface magnetic field is highly non-dipolar(for review see Geppert, Rheinhardt, 
& Gil 2003; Urpin & Gil 2003). The total surface field, S^, may be written as 

= bBa, (3) 

where 6 is a constant of proportionality which relates the dipolar field component, B^, to 
the total surface field strength. For VG models that include the effects of non-dipolar fields 
(i.e. 6 7^ 1), it will be useful to estimate b in terms of P and P. Two assumptions will be 
made in order to do this. First, it will be assumed that the observationally inferred surface 
magnetic field is actually a measure of the surface dipole component of the field. This is a 
reasonable assumption since the power radiated by the dipolar magnetic field dominates the 
other multi-pole components. The surface dipolar magnetic field strength, Bd, is estimated 

by 

Ba = 6.4 X 10^^(PP)°-^G = 2 x 10^^(PP_i5)°-^G, (4) 
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where P is the pulsar period and P is the period derivative with respect to time (Shapiro 
& Teukolsky 1983). The second assumption will be that the total surface field, S^, will be 
approximately equal to lO^^G for all pulsars, independent of P and P(Gil & Melikidze 2002, 
and references therein). With these two assumption, h may be written as 

6 = 5Si3(PP-i5)-°', (5) 

where Bi^ — S^/IO^^ G is of order unity. 

The size of the polar cap is determined by those magnetic field lines which connect to 
the interstellar medium magnetic field. For all the VG models discussed below, the polar 
cap radius will be estimated using 

Tp = b-^-^Rl^lO'^p-^-^cm, (6) 

where Rq is the neutron star radius in units of 10^ cm. This is the classical dipolar field polar 
cap radius modified by a factor of b~'^^^ due to magnetic flux conservation over the tube of 
open magnetic field lines. 

Four VG models will be discussed in the next four subsections. The models are classified 
by the structure of the surface magnetic field and the physics of the electron-positron pair 
plasma production process. The prefix "CR" stands for curvature radiation induced pair 
production (Erber 1966, RS75) while the prefix "ICS" stands for resonant inverse Compton 
scattering induced pair production(Zhang et al. 1997, and references therein). The suffix 
"VG" is used for models that assume 6 = 1 (i.e. a pure dipolar field) while the suffix "NTVG" 
is for "near threshold vacuum gap" models that assume b is given by Equation 5. The results 
of the next four subsections are summarized in Table 1 which gives a and a for each model. 



Table 1. Summary of Vacuum Gap models 



Model 


Complexity Parameter (a) 


a 


References 


CR-VG 
CR-NTVG 

ICS-VG 
ICS-NTVG 


3.07J-2/^R^/2s-//"p-w/28pl/t 

I.m6«-5^P^/'P0-I4p0.79 
4.57eg''"i?^/'B035p-0-39p0f| 


-2.25 
-2.71 

0.18 
-1.56 


1,2 
3,4 

3,5 
3,4 



References. — (1) Ruderman & Sutherland (1975), (2) Gil & Sendyk (2000), (3) Gil & Mitra (2001), (4) Gil & Melikidze 
(2002), (5) Zhang et al. (1997) 
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2.1. CR-VG model 

The model of RS75 is the prototype of a curvature radiation induced vacuum gap. The 
gap height is determined by the condition h = Iph, where Iph is the mean free path for a 
photon to coUide with the strong magnetic field and form an electron-positron pair. This 
model is valid for a dipolar magnetic field with < x 10^^ G. The gap height in this 
model is /i = (3.5 x lQ^)n^J^ P^/'^ PzlC cm (RS75, Gil & Mehkidze 2002), where 7^6 ^ 1 
is the radius of curvature of field lines in units of neutron star radius it! = 10^ cm. The 
complexity parameter, a, takes the form a = 2.97?.g ^^^i?g'^^P"^/^''P^(5 and a — —2.25. This 
value of a is currently in the range allowed by the analysis of JG03. 

2.2. CR-NTVG model 

The CR-VG model described above has a fundamental problem: the binding energy is 
much too small to prevent thermionic emission from the polar cap surface (e.g. Abrahams & 
Shapiro 1991; Usov & Melrose 1995, 1996). Gil & Mitra (2001) argued that this problem can 
be solved by assuming that pulsars have an extremely strong, non-dipolar surface magnetic 
field whose magnitude,!?^, is close to about lO^^G, independent of P and P. In such a 
strong magnetic field {Bs > 4:A x 10^^ G (e.g. Usov & Melrose 1995)), the curvature radiated 
photons will pair produce near the kinematic threshold (e.g. Daugherty & Harding 1983). 
Hence, this type of model is called a curvature radiation induced near threshold vacuum gap 
(CR-NTVG) model. The gap height in this model is /i = (3 x 10^)111^'^ b'^/^ P^/^'^PZ^i^'^ cm 
(Gil & Melikidze 2002, GM02 hereafter). Using Equations 1, 5, and 6 one obtains a — 
3.onf/'RTB;i/''p-'^/''p'_{t. This gives a — —2.71, amazingly close to the maximum 
significant value obtained by JG03. 

2.3. ICS-VG model 

In this model, the electron-positron pairs are produced by inverse Compton scattered 
seed photons in a relatively low strength dipolar magnetic field (Zhang et al. 1997, and 
references therein). The gap height in this model is h — (8.8 x iQ^^^j^o-br p-o.64:p-o.79 
(Zhang, Harding & Muslimov 2000; Gil & Mitra 2001). Thus, the complexity parameter is 
a = l.lT^g °'^''_Rg-^P°-^^P°/| and a = 0.18. This value is well outside the range of a values 
supported by the analysis of JG03. 
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2.4. ICS-NTVG model 

This is a near-threshold version of the ICS-VG model, valid for Bg > 4:A x 10^^ G. The 
gap height in this model is /i = (5 x 10)^7^g^'^6"^P"°-^^Pr°5^ cm (GM02). Using equations 
1 , 5, and 6 one obtains a = A.^Uq^-^'^ R^q^ 3^ p-^-^^ P^^l and a = -1.56. This value of a is 
also outside the current range of supported values. 

3. SUMMARY AND DISCUSSION 

By analyzing the pulse-to-pulse intensity fluctuations of a set of 12 pulsars, JG03 found 
a significant correlation between the measured intensity modulation index (Eq. [2]) and 
a{a) — P°'P when a is within the range [-5,-2]. The most significant correlation was found 
when a — —2.7. JG03 concluded that this modulation index relationship (MIR) is consistent 
with the RS75 vacuum gap model which predicts a — —2.25. The work presented here 
calculates the expected values of a for four models of the vacuum gap including the RS75 
model. Assuming that future observations confirm the results of JG03, the inverse Compton 
scattering (ICS) driven gap models can be ruled out since they predict values of a outside 
the allowed range. Note that JG03 performed their analysis on "core" emission components 
only, hence the ICS models are only ruled out for such components. The curvature radiation 
driven gap models are consistent with the JG03 results. Moreover, we demonstrate that 
the near threshold vacuum gap model of GM02 reproduces the observed value of the most 
likely exponent a — —2.7. Since this model requires large amplitude, non-dipolar, surface 
magnetic fields, confirmation of a MIR with a — —2.7 could be evidence for non-dipolar 
surface fields of order lO^^G. Such fields may be generated by small-scale turbulent dynamo 
action just after the formation of the neutron star (Urpin & Gil 2003) or by Hall-effect 
processes occurring on the stellar surface (Geppert, Rheinhardt, & Gil 2003). 

Part of this research was performed at the Jet Propulsion Laboratory, California Insti- 
tute of Technology, under contract with the National Aeronautics and Space Administration. 
GJ acknowledges the support of the Polish State Committee for scientific research under 
Grant 2 P03D 008 19. We thank E. Gil for technical help. 
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